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We consider an original variational approach for building new models of quintessence interacting 
with dark or baryonic matter. The coupling is introduced at the Lagrangian level using a varia¬ 
tional formulation for relativistic fluids, where the interacting term generally depends on both the 
dynamical degrees of freedom of the theory and their spacetime derivatives. After deriving the field 
equations from the action, we consider applications in the context of cosmology. Two simple models 
are studied using dynamical system techniques showing the interesting phenomenology arising in 
this framework. We find that these models contain dark energy dominated late time attractors with 
early time matter dominated epochs and also obtain a possible dynamical crossing of the phantom 
barrier. The formulation and results presented here complete and expand the analysis exposed 
in the first part of this work, where only algebraic couplings, without spacetime derivatives, were 
considered. 


CONTENTS 


I. Introduction 


0 


II. Relativistic fluid interacting with a scalar field: derivative coupling 

A. Lagrangian formulation and field equations 

B. Conservation equations 

C. Derivation of the conservation equations 


2 

] 

5 


III. Cosmology 

A. Cosmological equations 

B. Cosmological dynamics 

1. Model A: 2D dynamical system 

2. Model B: 3D dynamical system 

C. A glimpse at the perturbations 

IV. Discussion and conclusion 
References 


G 

G 

7 

8 

1C 

13 


E3 


I. INTRODUCTION 

The aim of the present work is to expand and complete the analysis performed in pj, which hereafter will be referred 
to as Part I of this study. In Part I new models of interacting dark energy have been introduced starting from a 
variational approach. Defining a consistent variational set up for interacting dark energy models is an important issue 
not only for exploring new physical models in late time cosmology, but also for producing new theoretical mechanisms 
responsible for the phenomenology of such interactions. From the observational perspective a coupling between dark 
energy and dark matter might give rise to some effects which could be detected in forthcoming surveys. In fact, 
although non-interacting dark energy models can still accommodate the observations [2j, such interaction is mildly 
favoured by the data 0], and can even alleviate some tensions present in the comparison of different datasets [4[. 
For these reasons the theoretical development of interacting dark energy models should proceed to both drive future 
observations toward possible detectable features and to justify all the phenomenological models considered so far in 
the literature; see e.g. (BJ and references therein. 
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The particular advantage arising from a variational formulation consists in generating fully covariant equations of 
motion, which can then be applied to cosmology at both the background and the perturbation level. It is a well 
known problem to understand how to promote the phenomenological background equations of interacting dark energy 
to their covariant counterparts, and a satisfactory solution is still missing [g], In fact, it will always be possible to 
construct two covariant theories which are equivalent at the background level but will have different perturbations. 
This is due to the non-linearity of the Einstein field equations in general. Hence, there is a strong motivation to 
construct models at the level of the action. 

A variational approach, such as the one advanced in Part I, permits to solve this problem by automatically producing 
the desired field equations. In Part I a canonical scalar field describing dark energy, namely quintessence, has been 
coupled to a barotropic fluid at the Lagrangian level; see also Q for similar ideas. The Lagrangian description of a 
relativistic fluid that has been used in Part I, and which will also be employed here, is the one outlined by Brown in 
[S'). In what follows we will not review the features of this formulation but refer the reader to Part I. All details and 
the thermodynamics are thoroughly discussed in Brown’s original paper. 

The main issue developed and studied in this Part II, is the possibility of directly coupling the 4-velocity of 
a relativistic matter fluid to a scalar field <fi. The simplest way of mixing U M and <fi in an interacting Lagrangian 
is through the scalar U^d^cj) which inevitably involves the use of a spacetime derivative. In Part I we considered 
only algebraic couplings between the scalar field and the fluid’s degrees of freedom where no spacetime derivatives 
are allowed. This in particular implies that the coupling U^d^cj) falls beyond the interacting quintessence models 
introduced and investigated there. The scope of this Part II is exactly to complete the analysis of Part I by studying 
derivative couplings , i.e. coupling where spacetime derivatives are allowed. 

The paper is organized as follows. In Sec. [TT| we will introduce the Lagrangian of a scalar field interacting with a 
relativistic perfect fluid through a derivative coupling. Brown’s formulation will be employed to define and treat the 
coupling terms which will be assumed to depend on as few spacetime derivatives as possible in alignment with the 
effective field theory approach. In Sec. IIIII we will then focus on cosmological applications, finding first the relevant 
equations governing the evolution of the universe at large scales, and then analysing them with dynamical systems 
techniques. Finally in Sec. [TV]we will discuss the results obtained, highlighting the differences between the algebraic 
couplings of Part I and the derivative couplings of Part II. In Sec. [IV] we will also draw conclusions and speculate on 
possible future works. 

Notation and conventions: Unless otherwise specified we will assume standard general relativistic notation with 
the metric convention (—1,1,1,1) and Greek indices running from 0 to 3. Sometimes the comma notation for partial 
derivatives will be used: for example q Units where c = h = 1 will be employed together with the definitions 
k 2 = M" 2 = 8t tG. 

II. RELATIVISTIC FLUID INTERACTING WITH A SCALAR FIELD: DERIVATIVE COUPLING 

In this section we generalize the algebraic couplings considered in the first part of this work [l] by allowing also 
terms depending on the scalar field’s first derivative d^cf) to interact with the fluid’s degrees of freedom. In analogy 
to the effective field theory approach, we will only consider couplings with as few spacetime derivatives as possible. 
Moreover, we will not consider non-canonical couplings to the scalar field’s kinetic term as this is beyond the scope of 
the present work. However, it would be a possible further extension of the variational approach to also discuss such 
models. 


A. Lagrangian formulation and field equations 

The total action of our interacting dark energy system is [l| 

S = J d 4 X (Ugrav + -Cm + + Ant) j 

where the gravitational sector £ grav is given by the standard Einstein-Hilbert Lagrangian 

r _ V 9 d 

•C-grav — 2^2 Cl , 

with R being the curvature scalar with respect to the metric g^ and g its determinant. The Lagrangian of the scalar 
field is assumed to be in the canonical form 


( 1 ) 


( 2 ) 


■£-0 — V 5 




( 3 ) 
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where V is an arbitrary potential depending on (ft. Using Brown’s formulation [li, j§J the Lagrangian for the relativistic 
fluid Cm can be written as 


C M = - V-~9 P{ n i s ) + J M {<P,n + + Pau a ) , (4) 

where p is the energy density of the fluid prescribed as a function of n, the particle number density, and s, the entropy 

density per particle. The fields <p, 6 and (3a are all Lagrange multipliers with A taking the values 1,2,3 and a .a are 
the Lagrangian coordinates of the fluid. The vector-density particle number flux J^ is related to n as 

J» = V=gnU», \J\ = y/-g lu ,Ji*J v , n=-^L, (5) 

where t/^ is the fluid 4-velocity satisfiyng t/ M C/ M = — 1. For further details regarding Brown’s Lagrangian formalism 

for relativistic fluids we refer to jsj. 

Considering the dynamical degrees of freedom of the above Lagrangian, we note that the only two possible scalar 
coupling terms (up to total derivatives) that can be constructed containing only first order derivatives of the scalar 
field </> are dcj ) 2 = g^ v d^,(j)d v (j) and J^d^cj). The first possibility is nothing but the usual kinetic term for the scalar 
field, while the second one is a new coupling term that can only be considered within the fluid Lagrangian formalism 
treated in this work. In particular we will only study a linear coupling to J^d^cj) leaving higher order couplings aside 
at present. Note that within an effective field theory framework such higher order terms would naturally be neglected 
at first order. Therefore, the coupling term we will study in what follows is 

Ant = f(n, s, , (6) 

where / is an arbitrary function. This is the most general coupling term where only one spacetime derivative appears. 

The variations of © with respect to the Lagrange multipliers ip, 0, /3 a and Lagrangian coordinates of the fluid a A 
give the equations 


p : 

O 

II 
a. 
a. ~ 

(7) 

e : 


(8) 

/3a: 

P 

II 

o 

(9) 

a A : 

(JA3a), m = 0, 

(10) 


which are not modified by the coupling with the scalar field; see [lj,|sj:]. Eqs. © and © stand for the particle number 
conservation constraint and the entropy exchange constraint, respectively. These can be rewritten as 


V„(n £/**) = 0 and V^(s n U^) = 0 , (11) 

where is the covariant derivative with respect to g^. Eqs. © and (HOD determine the dynamics of the Lagrange 
multipliers but are not needed for our scopes and will not be considered further in what follows; see Q for more 
information about their physical meaning. The variations with respect to and s yield, respectively 


pUf, + + sd ' iM + f3 A a A ^ = ^n^-U^U v - fS^j q\ v , 


T = une M + 


d J-A 


where the chemical potential p and temperature T are defined by 


dp 1 dp 

/z = — and T =- — 
on n on 


( 12 ) 

(13) 


(14) 


Eqs. m and (1131) show how the four-velocity fluid decomposition and the temperature of the fluid depend also on 
the scalar field <f> due to the non-vanishing interaction. Similar equations are found also for the algebraic couplings 
considered in Part I. They modify the thermodynamics properties of the fluid reflecting the fact that an interaction 
with the scalar fields is now at play. 

Variation of © with respect to g^ give the following Einstein field equations 






(15) 
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where 


— p Q[iv T (p T P) U^Uv , 


fils 


rr i(int) 

[IV 


d^d v ij) - g^ 


df 

-n 2 -±U x d x <l>{gp V + UJJ v ) . 




with the fluid pressure defined as 


dp 

p = n— - p. 

on 


(16) 

(17) 

(18) 


(19) 


Interestingly, now the energy-momentum tensor of the interaction is orthogonal to the fluid flow, = 0 and it 

can be rewritten in a perfect fluid form 


— Pint + (pint + Pint) , (20) 

with vanishing energy density and pressure given by 

df 

Pint = 0 and Pint = -n 2 -^U x d X (t> ■ (21) 

The scalar field equation, obtained from the variation of 0 with respect to </>, is 

a<j>-V' + n 2 ?J-S7u.U» =0. (22) 

dn 


Note that whenever / does not depend on the particle number density n the interaction term does not give any 
contribution to the equations of motion since T^ t ' > vanishes and Eq. (l22l) reduces to the uncoupled Klein-Gordon 
equation. This is expected since in this case conditions 0 and 0 would constrain the interacting term 0 to become 
a total derivative and thus a boundary term. Note also that the interaction contributes always in the form n 2 df /dn 
which does not depend on n if one chooses / oc l/n. This case is particularly interesting since such a coupling would 
correspond to A nt = y/^W^U^d^fj) with W((f>) being a general function of (j) (ignoring the dependence on s). Its 
cosmological consequences will be studied in detail in Sec. IIIII 

Intuitively, this result is not surprising. Our coupling is based on the idea of allowing an exchange of energy between 
the scalar field and the matter via the particle number density n of that matter. Removing this dependence from the 
interaction Lagrangian corresponds to closing this channel of energy transfer, which, given that the fluid 4-velocity is 
conserved, is precisely why © becomes a boundary term in that case. 

The original gravitational field equations that have been derived in this section, namely Eqs. m and (122p . represent 
a completely new way to couple a relativistic perfect fluid to a scalar field. Even when compared to the interacting 
equations obtained in [lj , one can realize that these equations not only take into account the coupling between 0, n 
and s, but heavily depend also on the fluid 4-velocity. These equations can be used to build new models of interacting 
quintessence (and coupled inflation) possibly leading to unexplored phenomenological and theoretical issues. 


B. Conservation equations 

In the following we will briefly discuss how our formalism compares with the standard approach of introducing 
couplings at the level of the field equations (see [§] for a general discussion on this issue). One can rewrite the coupled 
field equations above in the more familiar form 

GW = fay + = Q u and = -Q „, (23) 

which is commonly used to couple two matter components in general relativity. For this purpose we define a new 
energy-momentum tensor for the fluid as T M!/ = + T‘^ t ' > with the energy density p = p + p - m t and pressure 

p = p T pint • The Einstein field equations (1151) then become 

G^ = k 2 fav+T^) , 


(24) 
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resembling the first of Eqs. (l23l) . Using the Klein-Gordon equation (l22l) and the energy-momentum tensor urn the 
conservation equation for the scalar field can be written as 

V'TW = -n 2 ^V x U x V v <f> = Q v , (25) 

which shows that the scalar field is not conserved due to the interaction with the fluid. 

To prove that also the fluid energy-momentum is not conserved, it is better to split its conservation equation into 
the parallel and perpendicular components to the fluid flow: 

- U U U X V^ X , (26) 


where is defined by 


hfiv — g^v + . (27) 

Then using Eqs. m-m and Eq. (fl2ll one can show that (see the next subsection) 

h^V x f Xv = h"n 2 ^y x U x V„<f> and = !Tn 2 ^V A f/ A V^, (28) 

which inserted back into Eq. (l26l) gives 

V‘% v = n 2 ^-V x U x d v <j> = -Q v . (29) 

This shows that the interacting approach developed in the present paper can be mapped back into the standard 
relativistic coupling between two matter components as given by Eqs. (1231) with an exchange vector Q M defined in 
Eqs. (l25l) and (l29l) . Note however that the pressure of the coupled fluid p now differs from the uncoupled matter 
pressure p in as much as it depends also on the scalar field through the interacting term. For example, for a pressure¬ 
less fluid, where p = 0, in general p ^ 0. Note also that for the derivative coupling considered here, the energy density 
of the coupled and uncoupled matter fluids coincide, namely p = p since pint = 0. 

For the covariant conservations of the “bare” one instead finds 


= 0 and 


h^\/ x T vX = 2 nU x V [x 


df tt rrv 




(30) 


where square brackets denotes anti-symmetrisation. This implies in particular that the “bare” matter energy- 
momentum tensor is not conserved in the presence of the coupling with the scalar field 

V m T^^0. (31) 

This is not unexpected since the Einstein field equations imply the conservation of the total energy-momentum tensor 
only. Any additional conservation equations need to be imposed separately. 


C. Derivation of the conservation equations 


The derivations of Eqs. (E51) and (13(71) are slightly involved and we show some details here. These calculations are 
similar to the ones appearing in the Appendix of [lj. Recalling that = 0 and using Eqs. CCD- we have 


U X V^ = t/ A V^ ( pg^ x + (p + P )U fl U x - n 2 ^-U a d a (t>h^ x 


..2 rrao ATT^\ 


= -U^p - (p + p)V^ - n ^ U a d a <t>U A a 
= ~U^n - - n 2 ^-U a d a (j)U x S7 fi (U X U^) 


dn 


dn 


dn 


dp^ u , rr , . 2 d f, 


= {U » n) + n -^ ua V<*^ u » 


(32) 

(33) 

(34) 

(35) 


( 36 ) 
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in agreement with the first of Eqs. (|28|) . Note that if we consider instead of in the passages above, we obtain 


U X W^T^ = 0, 

proving the first of Eqs. m- For the second of Eqs. C51) we have 

h v ^ x T vX = h v y S7 x T vX - ^V A (n 2 ^U a d a <t>h x 

For the first term we can use the relation (see the appendix of 0 ) 

K v V x T vX = 2nU x V [x {iiU tA ), 

which gives 

7Arr> . _ o„rrAv7_ ( ,.tt .\ i.i'wA ( „2^/ i 


h^ x T„ x = 2nt/ A V [A (/R/ M] ) - h^ x ( n 2 -±U a d a <t>h x 


(37) 


(38) 


(39) 


(40) 


Then making use of Eq. m we obtain 


/£V A f vA = -2 nU x [V [a V m] ^ + 8 V [a V m] 0 + V [A {p A V^a A )} + 2 nU x \7 [x 


- fS *) v "^ 


2 <9/, 


-KX7 x n -^-U a d a (f)h x 


= nf7 A V A 


dn 


n % U » UV -f S »)Vv<i 


(41) 


- nU x V L 


n ^ UxUV ~f 5x ) V ^ 


2 df. 


-KS7 X [n 2 -^U a d a (f>h Xv 


(42) 


= n 2 |^V A 17 A (V^ + U V U^A) 

= K ( n2 ^ xUX ) V ^’ 


(43) 

(44) 


where to go from line (14T1) to (l42|) we used the fact that covariant derivatives commute on any scalar and we applied 
Eqs. CD-©; while the passage from line (H21) to C3l) requires long but standard algebraic manipulations together 
with the use of Eqs. CD- This proves the second of Eqs. (12811 . Note that if the same calculations m gd are applied 
to T^ v one immediately finds 


h"V x T vX = 2 nU x V [x 


df n TT» 

n dT u ^ u 




in agreement with the second of Eqs. 


(45) 


III. COSMOLOGY 

In this section we study cosmological applications of the theoretical framework outlined in the previous section. We 
will first derive the background equations and then rewrite them in the form of a dynamical system. Subsequently, 
employing dynamical systems techniques, we will study the universe evolution given by two particular interacting 
models defined within our approach. 


A. Cosmological equations 

We will now assume a flat Friedmann-Robertson-Walker (FRW) metric 

ds 2 = — dt 2 + a(t) 2 [dx 2 + dy 2 + dz 2 ) , (46) 

where a(t) is the cosmological scale factor and (i, x, y, z ) are comoving Cartesian coordinates. We will also assume 
that all the dynamical quantities are homogeneous, i.e. they depend only on the cosmological time t. In particular we 
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will have that (f >, p, n , s will be functions of t only. Moreover taking into account comoving coordinates the perfect 
fluid 4-velocity becomes simply U 11 = (—1, 0, 0,0). 

We can then derive the following cosmological equations from the Einstein equations (fl5ll 

3 J^ = p+ ^ + V, (47) 

-^(2H + 3 ^)=p + ^-V-n^, (48) 

where H = a/a is the Hubble rate and the over-dot denotes differentiation with respect to t. The scalar field equation 
is instead 


d> + 3H<fi+V' -n 2 ^J-3H = 0, (49) 

on 

and is obtained from (1^1) . We will work with the standard exponential potential V = Vo exp(—A/t</>), leaving the 
analysis for different potentials to future work. We can immediately notice that the Friedmann equation (l47l) is not 
modified by the interacting term. This happens because the time-time component of (1181) vanishes for the background 

FRW metric, or equivalently because is orthogonal to the fluid flow as pointed out for Eq. ED- 
In a cosmological framework we will always have that Eqs. ED give the constraints 

h + 3Hn = 0 and s = 0 . (50) 

These equations tells us that the entropy density per particle is conserved through the universe’s evolution, while the 
particle density decays according to 


n oc a 3 , (51) 

which is expected from geometrical considerations. The two dynamical quantities ultimately appearing in Eqs. ED 
(ED are thus <j>(t) and a[t) only. 


B. Cosmological dynamics 


In this section we will examine the dynamics of a universe described by the cosmological equations (ED ED using 
dynamical systems techniques. We begin by introducing the canonical dimensionless variables [HEU 




X = 


K(f) 


V3H ’ ~ V&H ’ V V3H ' 

which give us the Friedmann constraint 

1 2 , 2,2 
1 = a + x +y . 

Using these we can rewrite the cosmological equations as the dynamical system 

x' = —^(3a;((iu — l)x 2 + (w + 1 )y 2 + 1 — w) — \/6(A(x 2 — 1) + A y 2 )) , 
V = -h(3((w - l)x 2 + (w + 1 ){y 2 - 1)) + v^A - H)). 




(52) 

(53) 

(54) 

(55) 


The quantity A is defined by 


A K 2 d f 

a = - h " 3S 


(56) 


Hence the dynamical system is not closed until one specifies the form of the function /. For particular choices of 
/ we can rewrite A as a function of the two dimensionless variables x and y. These special models do not increase 
the dimensions of the dynamical system which remains 2-dimensional. In general, however, if A cannot be written 
in terms of x and y only, the dimension of the phase space will increase. This does not hinder the use of dynamical 
system techniques, however, visualising the phase space might become more involved. An example is a coupling of 






the form / = —/o exp( 7 K</)/n with 7 a dimensionless constant. In this case A = —k/o exp( 7 «(/)/.ff which cannot be 
expressed in terms of x and y only. 

In the following we will consider one 2 dimensional and one 3 dimensional model, with two choices for / as displayed 
in Tab.lH where a , 7 and £ are dimensionless constants, with 7 algebraically related to £ by 


7 (3 - a )( w + 1) - 1 


(57) 



/ 

A 

Model A 

Model B 

-yp 1 / 2 - a V a /V3n 

(Ho 

htn 

ty 2a ( 1 

- X 2 - y 2 ) 1/2 ~ a 

t H 0 

^ H 


TABLE I. Choices of the interacting function for the dynamical systems analysis. 

These models are studied using the standard approach of dynamical systems analysis. The critical points of the 
respective models will be found, their corresponding eigenvalues will be computed, and these results will be interpreted 
within the context of cosmology. When possible the phase space of the dynamical system will be drawn together with 
few trajectories numerically computed for particular choices of the model parameters. 


1. Model A: 2D dynamical system 


In this subsection we analyse model A of the dynamical system (l5H) - (l55l) . This corresponds to choosing the function 
/ to be 


/ °c 


p l/2-aya 

•\/3 n 


In this case the dynamical system remains two dimensional, as A can be written in terms of x and y as 

A = (y 2a (l — x 2 — y 2 ) 1/2- “ 


(58) 


(59) 


with £ a constant parameter. The Friedmann equations can be rewritten to give the acceleration equation 

= tH 1 + w) + (w- l)x 2 + (1 + w)y 2 - 4^(1 - X 2 - y 2 )i~ a ] (60) 

which means we can define the effective equation of state (EoS) 


w eS = x 2 - y 2 + iu(l - x 2 - y 2 ) + \J^(xy 2a (l - x 2 - y 2 )* “. (61) 

Critical points correspond to solving the system 

x' = 0, y = 0. (62) 

For general values of a finding critical points is difficult, and the system becomes singular unless a lies in the range 
0 < a < 1/2, so in the following we will examine some specific choices of a lying in this range. The simplest such 
choice is a = 1/2. This corresponds to choosing / = t;VV /(3n), meaning A = £y. The background dynamics 
of this model corresponds to the one of a particular fc-essence scalar field analysed in EH- The phenomenology 
at cosmological distances is quite interesting and includes late time phantom dominated solutions with dynamical 
crossing of the phantom barrier. Moreover while the model considered in El always exhibits instabilities at the level 
of cosmological perturbations, the one constructed with the formalism of this work could be stable since the dynamics 
at the level of perturbations will be completely different. A more in depth discussion on such model, namely for the 
choice / = (;VV /(3 n), will be given in Sec. Hvl 

In the remainder of this section we will examine the case a = 0. We set w = 0 for simplicity, since other values of 
w are not relevant for dark matter interacting models. The critical points of the dynamical system are displayed in 

Tab. [Hi 

Properties of these critical points, including existence and stability can be found in Table IIIII The system has 
potentially up to four critical points depending on the values of the parameters A and £: 
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Point 

X 

y 

A± 

±1 

0 

C 

A 

%/6 


D 

V2£ 

\/3+2£2 

0 


TABLE II. Critical points of Model A with w = 0 and a = 0. 


• Point A±. These two points exist for all A and £. They are the standard solutions dominated by the scalar 
field kinetic energy, with effective equation of state of a stiff fluid w e g = 1. These points are either unstable or 
saddle points depending on whether the absolute value of A is less than y/6. 

• Point C. This point corresponds to a universe completely dominated by a scalar field. It exists only for A 2 < 6. 
It is stable when £ 2 < 2(3 — A 2 ), and a saddle node otherwise. The point describes an accelerating universe 
when the scalar potential is sufficiently flat, requiring A 2 < 2. 


Point D. In the limit £ —> 0 this point reduces to the origin. The energy density of this point is dominated 
both by the matter and the kinetic energy of the scalar field, with no scalar field potential energy density 

The effective 


contribution. It is the stable late time attractor only for negative £ satisfying A£ < — | y 1 + ^ 

EoS is w e ff = 0, the same as a matter dominated solution, however it should be noted that for general w it is 
not a scaling solution, generally w e s = 3w/(3 + 2£ 2 ). 


Point 

Existence 

Weff 

Acceleration 

Stability 

A- 

VA, £ 

1 

No 

Unstable node: A > — \/6 

Saddle node: otherwise 

A+ 

VA, £ 

1 

No 

Unstable node: A < a/6 

Saddle node: otherwise 

C 

A 2 < 6 

A 2 —3 

3 

A 2 < 2 

Stable node: £ 2 < 2(3 — A 2 ) 

Saddle point: £ 2 > 2(3 — A 2 ) 

D 

VA, £ 

0 

No 

Stable node: A£ < — § \J 1 + 
Saddle node: otherwise 


TABLE III. Stability of critical points of Model A with w = 0 and a = 0 

The phase space for this two dimensional model is simply the upper half unit disc. We show the phase space 
diagrams with example trajectories in Fig. [Q Fig. [2] and Fig. [3] The region of acceleration is indicated by the grey 
region. The shape of this region of acceleration is dependent only on £, and is independent of the parameter A. In 
Fig. [Q we make the parameter choice A = 4 and £ = —1. Trajectories begin at the stiff matter point A_. Many 
trajectories then undergo a short transient accelerating phase, before then decelerating, with some trajectories being 
drawn towards the saddle point A_|_. Finally the late time attractor D is reached, which has a matter equation of 
state. 

In Fig. [2] we choose A = 2 and £ = +1. Trajectories now start at either of the stiff matter points A + or A_. 
Trajectories are then drawn towards the saddle matter dominated point D. Trajectories are then attracted upwards, 
passing through the region of acceleration. They are then finally attracted towards the global attractor C, which is 
not accelerating for this choice of parameter values. 

In Fig. [3] the parameter values A = 1 and £ = —3/2 are chosen. Again trajectories begin at the stiff matter points 
A + and A_ before passing through the matter dominated saddle point D. Trajectories then enter the region of 
acceleration and end at the late time global attractor C, which in this case is accelerating as the scalar potential is 
sufficiently flat. 
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These models are able to accurately describe late time universe phenomenology, there are many parameter choices 
which result in a global accelerating attractor. They can also describe universes undergoing transient periods of 
acceleration. However as is typical in these models it breaks down at early times, as it begins at the stiff matter point 
A-, whose effective EoS w e g = 1 is not physically viable. 

There are a few key differences with these models if compared to the canonical scalar field case, which we recover 
by taking the limit £ —> 0. There is no longer a scaling solution in this model, however the origin O is deformed into 
the point D which behaves as though the universe is matter dominated. Moreover, unlike the origin O this point D 
can be stable for a variety of parameter choices. 



-1.0 -0 5 0.0 0.5 1.0 

x 


FIG. 1. Phase space for the dynamical system when £ = —1 and A = 4. The shaded region indicates where acceleration is 
present. 



FIG. 2. Phase space for the dynamical system when £ = 1 and A = 2. The shaded region indicates where acceleration is 
present. 


2. Model B: 3D dynamical system 

In this section we consider model B of the dynamical system (15411 - (155)) where we take the coupling function / to 
simply be proportional to 1 /n, as outlined in Tab. Q] In particular we take 


/ = 


£# 0 


(63) 


where £ is a constant and Hq is the Hubble parameter at an arbitrary fixed time. The function A is now simply 


A = £— 
^ H 


(64) 
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FIG. 3. Phase space for the dynamical system when £ = —1.5 and Awl. The shaded region indicates where acceleration is 
present. 


The resultin g dy namical system is no longer two dimensional, in which case we need to introduce a third variable, see 
for instance [Tj]. We take this to be 


H a 


Hn + H 


(65) 


which is compact and lies in the range 0 < z < 1. The resulting dynamical system is given by 


(3z(z — 1)(1 — tu + (in — l)x 2 + (1 + w)y 2 ) + \/6(£t:(l — x 2 ) — A (z — 1 )y 2 )) 


V = 


2(*-l) 

V 

'2(*-l) 


(3(2 — 1)((1 + w)(y 2 — 1) + (w — l)x 2 ) + V& x((z — 1)A + z£)) 


z' = - 2 ( 3(2 — 1)((1 + w)(y 2 — 1) + (w — l)x 2 ) + VGz^x). 


( 66 ) 

(67) 

( 68 ) 


The Friedman equations can again be rearranged to give the acceleration equation in terms of these new variables 
x, y and 2 . This time for the acceleration equation we find 


H 3 r \ 2 /-t \ 2 2F ^X i 

Jp = 2 H 1 + w ) - (! - +(1 + w)y - 

Hence in this model w e g is given by 


(69) 


W e ff = W + (1 — w)x 2 


(w + l )y 2 + 


2£ zx 
^6 1- z' 


(70) 


The variable 2 was chosen so that the phase space of the system is compact, this time the phase space will be a 
semi-cylinder of unit height. The dynamical system has a singularity on the plane 2=1, which corresponds to the 
Hubble parameter H —> 0, i.e. when a —> 00 , usually happening at t —> 00 . However, none of the critical points of 
the system lie on this plane. We display the critical points for the system in Tab IIVI 

The stability of the critical points is shown in Tab. [V] The system has at most seven critical points at any one 
time. The points O, A±, B and C lie on the 2 = 0 plane and have exactly the same coordinates as the canonical 
scalar field [10(. Moreover, the existence and acceleration properties of these points remain the same as the canonical 
case. This happens exactly because the contribution of the interaction vanishes for 2 = 0, as can be realised looking 
at Eqs. (® G0- However, none of these points can now be stable, the 2 = 0 plane is unstable in general and 
trajectories starting on this plane will not stay there. There are three new critical points in this model, given by: 

• Point D This point exists only for £ > 0. At this point the energy density of the universe is dominated by the 
scalar field potential energy, with zero matter and scalar field kinetic energy contributions. When this point 
exists, it is always the late time global attractor. Moreover it has an effective EoS w e g = — 1, hence acceleration 
is present and the behaviour of the cosmological constant is mimicked. 


• Point E_ This point only exists for £ > 0. The energy density of the universe at this point is dominated by 
the kinetic energy of the scalar field, with no potential or matter contributions. When this point exists it is the 
unique late time attractor. It also has the effective EoS w e g = —1. 
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Point 

X 

y 

Z 

O 

0 

0 

0 

A± 

±i 

0 

0 

B 

/?(1+™) 

Hi \A 1 + U ')( 1-U ') 

0 

V 2 A 

v 2 A 

C 

A 

a/6 

A-t 

0 

D 

0 

1 

A 

A+£ 

E± 

±i 

0 

V6 


TABLE IV. Critical points of Model B 


• Point E .|_ This point is similar to i?_, except it only exists for £ < 0. It is also dominated by the scalar field 
kinetic energy, and has effective EoS w e g = — 1. However unlike £L, this point is always a saddle point and is 
unstable in general. 


Point 

Existence 

Weff 

Acceleration 

Stability 

O 

VA, £ 

w 

No 

Saddle node 

A- 

VA, £ 

1 

No 

Unstable node: A > — %/6 

Saddle node: otherwise 

A+ 

VA, £ 

1 

No 

Unstable node: A < y/E 

Saddle node: otherwise 

B 

A 2 > 3{w + 1) 

w 

No 

Saddle point 

C 

A 2 < 6 

A 2 —3 

3 

A 2 < 2 

Saddle point 

D 

£>0 

-l 

Yes 

Stable node: A 2 <3/2 

Stable spiral: A 2 >3/2 

E+ 

£ < o 

-l 

Yes 

Stable node 

E- 

£ > o 

-l 

Yes 

Saddle point 


TABLE V. Stability of critical points of Model B 

The global dynamics of these models are particularly simple to analyse. The late time global attractor is either D, 
in the case £ > 0, or E+, in the case £ < 0. Moreover both of these late time global attractors have effective EoS 
w e g = — 1. Trajectories begin on the z = 0 plane, and while they remain on this plane the trajectories behave as 
in the case of the canonical scalar field. This plane is generically unstable and the trajectories eventually leave this 
plane and arrive at the late time attractor. 

Phase spaces for a couple of parameter choices with example trajectories are plotted in Fig. [4] and Fig. [5] In 
Fig. 21 the parameter values A = 2, £ = — 1 are chosen. Trajectories start at the stiff matter points A±, before being 
attracted towards one of the three saddles O, B or C on the z = 0 plane. All trajectories then end at the late time 
accelerating global attractor E + . There is a particular (heteroclinic) trajectory of interest that passes through the 
matter dominated point O before ending at the accelerating point E+. Ignoring the theoretical issues at early times, 
this trajectory would mirror the dynamics of a universe with no scalar field potential but still with a cosmological 
constant late time behaviour. 

In Fig.[5]the parameter values are set to be A = 2 and £ = +1. Trajectories start at the stiff matter points A±, and 
are then attracted towards the saddle points on the plane 0 = 0. The trajectories then leave the plane, with some 
attracted towards the saddle point E_. All trajectories end up at the late time accelerating attractor D. 
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X -10 



0.5 


FIG. 4. Phase space of model B of the dynamical system with values w = 0, A = 2, £ = —1. 


These models are able to accurately describe late time universe phenomenology. All trajectories end at accelerating 
critical points with effective equation of state w e g = — 1. In this scenario the dynamics of the universe would mirror 
that of a universe with no scalar field but with a cosmological constant and thus these models could represent a 
solution to the cosmological constant problem. 


x -io 



0.5 


FIG. 5. Phase space of model B of the dynamical system with values w = 0, A = 2, £ = 1. 


C. A glimpse at the perturbations 

In this subsection we will briefly consider scalar perturbations of Model B in the linear approximation. Note that 
from the dynamical systems analysis of Sec. IIII B 21 it is clear that the late time attractor of this model is always 
represented by a de Sitter solution, which is well suited to describe dark energy domination. Moreover, this is the 
simplest scalar-fluid coupling where the gradient of the scalar field (f> appears. It will thus be interesting to see how 
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the cosmological perturbation equations are modified in presence of such a coupling. In what follows we will show 
the perturbed equations in Newtonian gauge and briefly discuss their main features. A complete investigation of the 
dynamics of cosmological perturbations is outside the scope of the present work, but will be delivered in forthcoming 
studies [l3j . 

We will study the behaviour of scalar perturbations assuming the following metric in Newtonian gauge using 
Cartesian coordinates 


ds 2 = —(1 + 24>)d/ 2 H - —— a - g (dx 2 + dy 2 + dz 2 ) , (71) 

[l + \k(x 2 + y 2 + z 2 )\ 

where k = —1, 0,1, and 4/ and $ are functions of all the coordinates. Since in both the explicit and implicit frames 
all the matter sources can be written as perfect fluids, no anisotropies appear in the scalar-fluid models considered 
here. In fact, the off-diagonal spatial field equations immediately imply that 4> = 4q which will be used to simplify 
the other perturbed equations. We will also give the equations directly in the Fourier space: V 2 i-» —q 2 , where V 2 is 
the Laplace operator. 

The covariant field equations to perturb are given by Eqs. (USD (EH)- As we mentioned, we will restrict the equations 
to the particular coupling (1631) . and, in order to simplify the notation, we will also redefine the coupling constant 
as £ = £Hq/k. Now, p, p and </> will refer to background quantities, while 4q 6(f), Sp and Sp will denote small 
perturbations. Furthermore, the perturbed 4-velocity of the fluid is 

5U fl = (-*,d i v), (72) 

with v being the scalar perturbation of the matter fluid velocity. The time component of this relation, with 4q follows 
from the constraint = — 1. 

At this point we are ready to state the equations for the cosmological perturbations. The 00-conrponent of the 
perturbed Einstein field equations then reads 

(+ 8np + 87tF — 6—^ 4/ + 3774/ + AirV'Scf) + 47 r5p + Aird)S(f> = 0 , (73) 

\ a 2 a z ) 

while the 07-components are 


47t (£<f) + p + pj v — 47 j(t>8(j) + 774/ + 4> = 0 , 
and the diagonal ij-components become 


4/ + 477 'F + ^47 t£</> + 277 + 377 2 + 47n/> 2 - ^ 47 — 47T 1 5<j> + 4,TrV'5(f> — AnSp = 0 . 


The perturbation of the scalar field equation is given by 

„2 \ „2 


( 2 \ 2 

^ + F"J 5<j>-£^v+ (377|+2W) (3| + 4^4- = 0. 

Finally we also provide the perturbations of the matter conservation equations: the time-component is 

q 2 

( p + p)—^v — 3H5p — 3H5p + 3 (p + p)'F — Sp = 0 , 


(74) 


(75) 


(76) 


(77) 


while the spatial components are 


£ + 3Hcj)\ - 377c 2 (p + p) v + v (^4> + P + pj + 3£HSc/) + £S(/) + (p + p)4/ + Sp = 0 , 


(78) 


where c 2 = dp/dp. 

We should firstly notice that the coupling does not appear in Eqs. GSD and (ED. This is a general feature of these 
derivative scalar-fluid models as will be shown in [l3j]. Mathematically this is due to the fact that the interacting 

energy-momentum tensor is orthogonal to the fluid’s 4-velocity, namely U ^7)1”' t ' 1 = 0, as discussed after Eq. (fTTH) . In 
fact, in a cosmological context this implies that the time-component of both, the Einstein and matter conservation 
equations, reduce to their uncoupled counterparts, as already noticed for the background equations. 
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The remaining perturbation equations are modified by terms containing £, which is the only coupling parameter 
appearing in these models and thus fully characterizing the interaction between dark energy and matter. In the limit 
£ —> 0 Eqs. (1731) (1751) reduce to the uncoupled correspondent perturbation equations. Note also how the interacting 
terms mix the gradient of the scalar field and the velocity of the fluid. 

This is expected from the derivative coupling used in the Lagrangian, which thus represents interesting new inter¬ 
actions even at the level of perturbations. In [131 ] it will be shown that the evolution of the perturbations governed by 
Eqs. (1731) (1751) are indistinguishable from the ACDM dynamics at sub-horizon scales (q H 2 a 2 ), although signatures 
of the interaction should arise at larger or non-linear scales. 


IV. DISCUSSION AND CONCLUSION 

The main motivation of this paper was the continuation of the approach outlined in Part I by taking into account 
derivative couplings between a scalar field and a barotropic fluid based on a variational formulation. In particular, we 
were interested in couplings linear in the first partial derivatives of the scalar field. The general expression for such a 
coupling can be written in the form f(n, s, (f^J^d^ and can only be treated with the methods previously developed 
in Part I. The presence of the term requires a separate treatment from that in Part I since the effects given 

by the appearance of the scalar field’s derivative will in general produce a different phenomenology. 

As an example we recall the fact that the energy-momentum (1181) . arising from the interacting contribution, is 
always orthogonal to the fluid 4-velocity, namely WT^' 1 = 0. Formally this happens because oc h which in 
turn is due to the fact that in the derivative interacting term © the metric tensor g appears within n only. As a 
consequence vanishes whenever f(n,s,(f>) does not depend on n, as one can easily realize looking at Eq. (1181) . 

The property U^T^' 1 = 0 is particularly relevant for cosmological applications since it implies that at the background 
level the Friedmann equation is never modified by the interaction between the fluid and the scalar field, see Eq. 63- 
This is in contrast to the algebraic couplings considered in Part I where this always happens. Note however that 
both the acceleration equation (1481) and scalar field equation (l49l) do get modified by the interaction. This implies 
that the background cosmological evolution can actually differ from the corresponding non-interacting one, though 
the Friedmann equation forces a matter-like evolution whenever the scalar field energy density becomes negligible if 
compared to the matter energy density, irrespectively of the strength of the interaction. The fact that = 0 

implies also that the perturbed Einstein equation (1731) as well as the continuity equation (1771) are not modified by the 
interacting term, although in all the other perturbed equations deviations due to the coupling appear. 

In Sec. uni we have studied some of the cosmological consequences obtained from these models, employing in par¬ 
ticular dynamical systems techniques. One would expect that a derivative coupling would complicate the subsequent 
dynamical systems formulation. However, it turns out that the derivative models studied here are somewhat easier to 
handle than models based on the algebraic coupling f(n,s,<j>) considered in Part I. In particular we have considered 
two types of derivative couplings, labelling them as Models A and B according to Tab. Q] The background cosmolog¬ 
ical evolution of Model A has been characterized by a 2D dynamical system, while the one of Model B required the 
analysis of a 3D dynamical system. For both models we found late time accelerating attractors, capable of describing 
the present dark energy dominated epoch, and possible intermediate phases of dark matter domination. 

For Model A the specific case corresponding to the parameter a = 0 has been studied in depth. We found that a 
dark matter to dark energy transition can be obtained for some values of the model parameters and is always described 
by the (heteroclinic) trajectories connecting Point D to Point C as depicted in Fig. [3] Note that, as mentioned before, 
Point D describes a scaling solution only when w = 0, i.e. when the fluid describes (dark) matter. For a different 
matter EoS one finds w e a ^ w, implying that modifications to the standard cosmological evolution might arise if 
radiation is added into the analysis. 

Another simple case of Model A is achieved by the choice a = 1/2. As mentioned above, the resulting background 
cosmological equations corresponding to this choice match the one arising from a particular fc-essence model studied 
in [lT|. Since the dynamical analysis of the a = 1/2 case has already been performed in detail in [Tlj] . it has not 
been considered in the present work. However we briefly recall here some of the features and results derived in [Til ] 
which arise as well in the background cosmological dynamics of Model A with a = 1/2. The phenomenology found 
in this case is highly rich with different cosmological behaviours that cannot be obtained with a non-interacting 
canonical scalar field. At late times it is possible to obtain phantom domination (w e g < —1) with dynamical crossing 
of the phantom barrier, while at early times super-stiff (w e g > 1) solutions can be attained. Theoretical models 
predicting a dynamical crossing of the phantom barrier are important to investigate in the eventuality that future 
surveys will detect (with statistical significance) a dark energy EoS within the phantom regime. Models of such kind 
obtained within the fc-essence framework, such as the one considered in [llj] . suffer from instabilities at the level of 
perturbations [Til ]. Nevertheless if the same phantom crossing evolution is not realized by a fc-essence model, as in 
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the case of Model A, then the cosmological perturbations might be stable. This is because, although the background 
equations coincide, one expects differences in the perturbation equations between Model A with a = 1/2 and the 
model considered in 0- These considerations enforce the motivation for analysing the dynamics of cosmological 
perturbations arising from the quintessence interacting models studied here and in Part I [13[. 

Considering Model B, it was found that a global accelerated attractor always appears in the phase space, irrespective 
of the values of the model parameters. This implies that for this particular coupling late time dark energy domination 
is always attained. Interestingly one of these attractors (Point E+), although mimicking a cosmological constant 
behaviour, is completely dominated by the scalar field’s kinetic energy, showing that within Model B late time 
acceleration can also be achieved without a self-interacting (and thus massless) scalar field. Moreover for some choices 
of the parameters a scaling solution attracting all the early time trajectories (the ones on the z = 0 plane) is present 
in the phase space; see Point B in Fig. [5] This saddle point forces the early universe to reach matter domination 
before switching to the late time accelerated behaviour. Thus for all physically possible initial conditions (for which 
z ~ 0, i.e. H Hq) a transition from dark energy to dark matter is attained at late times, solving in this manner 
any fine tuning issue. Although the presence of a scalar field at early times is expected to be strongly constrained by 
observations [151 ] , the results obtained within Model B are promising and should merit further consideration. For this 
reason in Sec. 1III Cl we presented and briefly discussed the scalar cosmological perturbations for such model, although 
a more detailed investigation will be presented in future works fl5 ]. 

To better expose the dark matter to dark energy transition at late times, we consider the behaviour of the effective 
equation of state for particular trajectories within the phase space of Model B and compare its qualitative properties. 
For this purpose, we look at the dynamical evolution of w e g for different choices of the parameters of Model B; see 
Fig. [G] We see that in this model we always find physical trajectories starting with the usual scalar field kinetic 
energy dominated epoch and then evolving through a matter dominated phase, after which they approach the late 
time dark energy dominated attractor. Note that in this case there is the possibility of dynamically crossing the 
phantom barrier, after which the effective equation of state approaches w e g = — 1 from below; see the left panel of 
Fig. |6] However, this behaviour is parameter dependent. For instance, when choosing £ = —1 the value w e g = —1 is 
approached from above. Note also that a non-monotonic transition from matter to dark energy domination can be 
obtained, as shown again in the left panel of Fig. [6] All this shows the interesting phenomenology present in Model B 
and suggests further investigations with other similar models. 

We might for example compare the coupling appearing in Eqs. (1251) and (1291) with the one usually considered in the 
dynamics of early-universe bubble nucleation arising from first-order phase transitions; see e.g. |l6j . In those models 
the coupling vector Q M between the scalar field and the surrounding matter fluid is provided by Q M oc U v d v (j)d^(f> and 
it is motivated by thermodynamical properties of electro-weak physics, such as the temperature dependence of the 
Higgs potential under renormalization at two or more loops. This coupling vector is different from the one arising 
in Eqs. (l25l) and (l29l) which assumes the form Q M oc V V U V d^cf) (and also retains a quite general dependence^] on the 
coupling function /). Unfortunately it seems not easy to obtain the coupling Q M oc U u d^cftd^cj) from the scalar-fluid 
variational approach considered in this work, even if higher-order operator such as are taken into account in 

the Lagrangian. Nevertheless, at a more phenomenological level, the coupling vector Q M oc V, y t/"<9 A1 </> can be used 
to describe new dissipative interactions between the expanding bubble and the surrounding fluid, which might give 


1 In EDI it will be shown that only models where / oc 1/n are generally free from instabilities at the cosmological perturbations level. 
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rise to a different dynamics in the early universe and possibly to new observational signatures of first-order phase 
transitions. These topics might well represent future applications of the variational formalism developed here. 

Finally we briefly comment on the appearance of a “fifth” force acting on the matter fluid due to the interaction 
with the scalar field. From Eq. (I4dl) (or equivalently from Eq. (1451) 1 it is possible to derive the geodesic equation for 
the fluid which reads 


dU» 

dr 


r Zpiruf 


w v 


p + p + Pint 


2 d f, 


d u p + <9„Pint - \U X W p 

on 


r, 


(79) 


where p- ln t is defined in Eq. m- Note that this fifth force is orthogonal to the fluid flow, / M t/ M = 0, in agreement 
with the relativistic definition of 4-force. The presence of a non-vanishing force on the right hand side of the geodesic 
equation, even for a dust fluid (p = 0), implies that in general the motion of matter particles will be non geodesic and 
that the equivalence principle will be violated. If the scalar field interacts with baryonic matter, then Solar System 
experiments set strong constraints on the magnitude of this fifth force E3, which must then be somehow negligible 
at small scales. However since it depends on local quantities such as the scalar field, the particle number density 
(or energy density) and the entropy density, a screening mechanism, similar to the well known chameleon mechanism 
[18l | . could be at play. The interesting and original feature of the non geodesic force (1791) is that it also depends on 
the fluid 4-velocity U M . This suggests the possibility of finding new screening mechanisms for the scalar field which 
hide its effects wherever the matter velocities are relativistically small, such as in the Solar System. Note that this 
cannot be achieved neither with the usual chameleon theories nor with the algebraic couplings of Part I, since in 
that cases the resulting fifth force does not depend on the fluid 4-velocity. Unfortunately the resulting analysis at 
small scales is complicated by this dependence on the matter velocities and cannot be easily performed following the 
original chameleon work, as it has been done for the algebraic couplings in Part I. The study of these new screening 
mechanisms depending on the matter fields velocities, as well as their phenomenological effects, falls well beyond the 
scope of the present paper and will be left for future considerations. 

In conclusion the variational approach outlined in this work is very powerful because it always allows us to arrive 
at covariant theories where the background and the perturbations can be studied consistently. It can also represent 
the starting point for more radical interactions between a scalar field and a matter fluid. For instance, one could 
attempt to identity the scalar field <j> with either of the two Lagrange multipliers ip or 9. In particular the first 
one would allow for fluid’s particle creation/annihilation, while the second one would correspond to the scalar held 
affecting the entropy of the system. Clearly, there is a large number of possible realisations whose consequences and 
applications could be studied in the future using different techniques, such as dynamical systems and perturbations. 
The investigation of cosmological perturbations constitutes the next logical step for further understanding the physics 
of these models and to better compare them against observations. 
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